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During inflation, primordial energy density fluctuations are created from ap- 
proximate de Sitter vacuum quantum fluctuations redshifted out of the hori- 
zon after which they are frozen as perturbations in the background curvature. 
Q_i! In this paper we demonstrate that there exists an intrinsic theoretical uncer- 

tainty in the inflationary predictions for the curvature perturbations due to 
the failure of the well known prescriptions to specify the vacuum uniquely. 
Specifically, we show that the two often used prescriptions for defining the ini- 
tial vacuum state - the Bunch-Davies prescription and the adiabatic vacuum 
prescription (even if the adiabaticity order to which the vacuum is specified 
is infinity) - fail to specify the vacuum uniquely in generic inflationary space- 
times in which the total duration of inflation is finite. This conclusion holds 
despite the absence of any trans-Planckian effects or effective field theory cut- 
off related effects. We quantify the uncertainty which is applicable to slow 
roll inflationary scenarios as well as for general FRW spacetimes and find that 
the uncertainty is generically small. This uncertainty should be treated as a 
minimal uncertainty that underlies all curvature perturbation calculations. 
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1 Introduction 



Recent measurements of the Cosmic Microwave Background (CMB) temperature fluctu- 
ations [1] give strong support for the picture that inflation gave rise to a flat spatial ge- 
ometry and the scale invariant energy density fluctuations on superhorizon scales needed 
for large scale structure formation. As a generic prediction of inflation, the primordial 
scale invariant energy density fluctuations on superhorizon scales can be calculated per- 
turbatively once the model of inflation is specified, including the prescription for the 
vacuum (i.e. Fock basis) of the inflaton field. 

To see what one means by specifying a vacuum in the canonical formalism, consider 
the simple case of a quantized scalar field (f> in an FRW-type universe with the metric of 
the form 

ds 2 = a 2 (r)(dr 2 -dx 2 ) (1) 

(conformal time coordinates) where we will restrict ourselves to flat spatial sections 
for simplicity. Without any nongravitational interactions, the field in the Heisenberg 
representation has an expansion 
r d 3 k 

0(x, t)=J (27r)3/2a(r) [a k h k (r)e ik - + al^K^] (2) 

where the a k is the annihilation operator which annihilates the vacuum and defines 
the Fock space. Because the creation and annihilation operators obey the commutator 
[ofcijajy = S^(ki — k 2 ), the h k s obey a normalization condition h k h' k * — h' k h* k = i to 
satisfy the canonical field commutators (henceforth, all primes or dots on functions of 
r refer to derivatives with respect to r as usual). Because the Heisenberg equation of 
motion for <fi will force h k (r) to obey a second order ordinary differential equation, it 
has two independent solutions, and hence there needs to be a boundary condition or 
prescription to specify h k thereby defining the vacuum state. 



Following the procedure in Minkowski space is one way of defining the vacuum (i.e. 
specifying hk)' define hk to be the positive frequency eigenvector of the time translational 
Killing vector. Because cosmologically interesting spacetimes do not have timelike Killing 
vectors, there is no a priori unique definition in specifying the vacuum even without any 
non-gravitational interactions. Indeed, even with a timelike Killing vector, as in de Sitter 
space, there can be some ambiguity [2, 3, 4, 5]. More generally, the simple reason for a 
lack of a good definition of vacuum state is just a problem of lacking asymptotically free 
states. The classic studies of quantum fields in curved spacetime (for good reviews, see 
for example [6, 7, 8, 9] and references therein) had focused primarily on systematically 
identifying the ambiguities and exploring various prescriptions that can remove them 
whenever possible. 

Two widely used prescriptions for defining vacuum are what is usually called the 
Bunch-Davies prescription [10] and the adiabatic vacuum prescription of Parker and 
Fulling [11, 12, 8]. The main point of this paper is to show that both prescriptions, even 
when the adiabatic order of the vacuum construction is infinity (i.e. the best one can 
do), do not specify the vacuum uniquely in most situations. This is true as well in one of 
the most important case of generic inflationary spacetimes 4 in which the total duration 
of inflation is finite, despite the absence of any trans-Planckian effects or effective field 
theory cutoff related effects [13, 14]. (For arguments against the trans-Planckian and 
cutoff effects, see [15, 16, 17, 18].) The simple reason is that both methods rely on 
an asymptotic definition of vacuum, leaving an ambiguity inherent in any asymptotic 
expansion. Note that the ambiguity that we are focusing on is also independent of 
the effects due to transition into inflation explored by Ref. [19, 20, 21]. In addition to 

showing the existence of an ambiguity, we estimate the uncertainty and its implication 

4 Hcrc, the "vacuum" of the inflationary era is referring to a no particle state of curvature 
perturbations. 
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for inflationary prediction of density perturbations. 

This work should complement the recent efforts [13, 14, 21, 22, 23] to uncover small 
quantum effects for the CMB measurements. After all, to see whether small effects can be 
measured, one must understand the inherent theoretical uncertainty in the calculations. 
Our uncertainty should be viewed as a minimal uncertainty that underlies all of these 
calculations, if one accepts either the adiabatic vacuum formalism or the Bunch-Davies 
vacuum formalism. 

In most generic situations, the uncertainty is extremely small. Indeed, in view of our 
work, it should be clear that the vacuum used in the work of Ref. [14] does not consti- 
tute the best adiabatic vacuum possible in a realistic inflationary scenario without dS 
invariance. From the adiabatic vacuum formalism, it should be considered an "excited" 
state. Even from a Hamiltonian minimization point of view of Ref. [16], one reaches 
a similar conclusion. To keep our presentation short as possible and to emphasize its 
independence from trans-Planckian issues, we will not discuss this point further in this 
paper. However, the reader should note that the existence of a cutoff does not change 
any of the results in this paper for the adiabatic vacuum. 

The order of presentation is as follows. We start off by reviewing the physical reasons 
for the ambiguities of a quantum vacuum in a cosmological spacetime. We then define 
and estimate the generic uncertainties of the "Bunch-Davies" and the adiabatic vacuum. 
The following section gives an example of how nonperturbative quantities can be calcu- 
lated (despite the usual inherent uncertainties) in an adiabatic vacuum in a very special 
limiting situation. In Sec. 5, we compute what can be seen as a slightly more precise 
uncertainty in the adiabatic vacuum during slow roll inflationary spacetimes. Finally, 
we summarize and conclude. 
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2 Ambiguities of vacua 

From a traditional particle physicist's point of view, vacuum can be defined as the state 
of no real particles. As has been well studied since 1960 's (see for example [6, 7, 8, 9] and 
references therein), this notion of vacuum is well known to be ambiguous: particles cannot 
always be unambiguously defined in the presence of a background field. For example, 
suppose one defines particles as the eigenstates of the momentum operator. These states 
as wavefunctions must then necessarily be spacetime translation eigenstates. However, if 
there is no spacetime translational symmetry of the background spacetime, there cannot 
be such an eigenstate. 

Even if the spacetime were flat Minkowski space, if one were to define particles 
empirically with an idealized detector, whether or not the detector registers particles 
depends in general on the motion of the detector. If the detectors were restricted to 
geodesies of the background spacetime, only in Minkowski spacetime, would all the 
geodesic detectors agree to no particle detection [8]. 

Because of this property of the Minkowski space, one may then try to argue that it is 
best to abandon the notion of a vacuous curved spacetime and treat it as a collection of 
gravitons in a Minkowski background. In this case, by definition the curved spacetime 
is not vacuous, although the state with background gravitons may be considered to be 
vacuous of some other field, say the inflaton field 0. Unfortunately, partly due to in- 
teractions with the background graviton fields, one is again forced back to having some 
ambiguities in defining the vacuum for 0. Hence, seen this way, we see the problem of 
the ambiguity of the vacuum, say with respect to the field 0, is not special to curved 
spacetime geometry, but to any situation in which there is a background field that in- 
teracts with the field 0. Just as in any other interacting quantum field theory, one may 
try to define free asymptotic states for and the graviton, and treat the interactions 
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perturbatively. This would work if there is an asymptotically flat region of spacetime. 
Unfortunately, the key difference in the curved background situations of interest to in- 
flationary cosmology is that there are no such asymptotically flat regions. 5 Hence, the 
problem of vacuum prescription can be rephrased as trying to define a vacuum state in 
the absence of free asymptotic states. 

Note that instead of relying on a particle definition of vacuum which is inherently 
nonlocal and observer dependent (the number of particles detected by any physical de- 
tector can be zero in one frame but not in another, even in Minkowski space), one can 
characterize the vacuum in terms of vacuum expectation value of the stress tensor, which 
has the advantage of being covariant (if it is zero in one frame, it is zero in all frames of 
reference). However, even here, the stress tensor vacuum expectation value is sensitive 
to the ambiguities in the boundary conditions and renormalization prescriptions of the 
correlation functions. These boundary condition sensitivities then reflect the ambiguities 
of the vacuum. 

3 The "Bunch-Davies" Vacuum and the Adiabatic 
Vacuum 

In this section, we would like to explain the two commonly used prescriptions that will 
be the focus of this paper. 

3.1 "Bunch-Davies" vacuum 

In the FRW cosmo logical context, the most well known and appealing prescription for 
identifying the vacuum is what is commonly called the Bunch-Davies prescription, which 

states that the positive frequency mode function (see Eq. (2)) should match asymp- 

5 There may be a way of obtaining a physically sensible answer by artificially turning on and off the 
gravitation, but we will not pursue this line of reasoning in this work. 
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totically to the Minkowski space prescription in the limit that the physical momentum 
k/a is much larger than the background geometry curvature scale (H = a{r)/a 2 ): 

lim h k ~ -^=e- ikT . (3) 
k/(aH)^oo y/2k 

In some cases, this uniquely specifies the vacuum. 

However,when this prescription is carried out at a fixed cosmological time (i.e. k — > oo 
with a(t) fixed), then the vacuum is not unique even in the simplest situations. Since 
inflation generally did not last an infinitely long time, the Bunch-Davies prescription 
cannot be applied in the asymptotic past limit, and the ambiguities associated with 
k — > oo are nonvanishing. This is independent of the existence of cutoffs or any trans- 
Planckian physics. 

For example, consider a massless scalar field minimally coupled to Einstein gravity 
in a patch of de Sitter space 

S = J (0 2 - (V0) 2 ) . (4) 

where a(r) = —1/(tH) > with l/H > being the dS radius. The resulting mode 
equation is 

h k (r)+wl(r)h k (T)=0, (5) 
where dots stand for derivatives with respect to the conformal time r and 



wl = e-* = k>-?- 2 . (6) 



The general positive frequency mode solution is 

„ e-' lkT ( iH \ „ e lkT ( iH \ 
V2k V (k/a) J V2k \ (k/a) J 

where A k and B k satisfies 

\A k \ 2 - \B k \ 2 = l (8) 



due to the normalization conditions given below Eq. (2). At any given time r, one 
can impose that the modes become Minkowskian as \(k/a)/H\ — > oo. Explicitly, this 
amounts to a functional matching of the form 



h k V2ke tkT -> 1 (9) 
for as \kr\ — > oo for all r < 0. Now, note that since r — > — oo in Eq. (7) gives 

h ^ A "7S + Bt Wk (10) 



for all > 0, this uniquely specifies A k — 1 and B k = for all A; > 0. 

However, consider the situation in which we make the restriction \r\ < |t | where To 
is the time at the "beginning" of inflation. Now, the limit \(k/a)/H\ — > oo can be taken 
only by taking k — > oo with a > a(r ). This means that we only know, for example, that 

B k - G(l/A; n ) (11) 

with n > as k — > oo. 6 Therefore, if we just impose /c — > oo with a bounded a(r), the 
Bunch-Davies prescription does not uniquely specify a vacuum. Applied to the case of 
Eq. (7), one must allow the ambiguity 

B k <0 (t^tt) ■ (12) 
\k/a(r )J 

A useful quantity to characterize the properties of the quantum perturbations of a mass- 
less scalar field during inflation is the power spectrum. For a generic quantity #(x, r), 
which can expanded in Fourier space as 

the power spectrum can be defined as 

<0Kflb|0) = 5 (3) (h - k 2 ) ^V g (k), (14) 



6 Note that exp (— ikr) does not have an asymptotic expansion in real k as k — > oo because of an 
essential singularity. 



where |0) is the vacuum quantum state of the system. This definition leads to the usual 
relation 

<(% 2 (x,t)|0> = / f V g (k). (15) 
If we compute the variance of the perturbations of the field 



<O|(0(x,r)) 2 |O> = / ^ |0 fe l 2 

dk k 3 



r dk k u 2 
J T W 1 kl 

/fnw, as) 



we may infer the power spectrum of the fluctuations of the scalar field to be 



k 3 



WO = 2^2 1^1- ( 17 ) 

Therefore, for a massless scalar field in de Sitter space, we obtain on superhorizon scales 
the power spectrum 

with = 1. From Eq. (18) we infer then that any ambiguity in the parameter B k 
implies an ambiguity in the power spectrum V<j, of the form 
SV^k) 



2\ReB k \~0(e~ (N »- Nk) ), (19) 



where N k denotes the number of e-foldings before the end of inflation when a given 
wavelength A = a/k leaves the horizon during the de Sitter stage and iVo denotes the 
number of e-foldings before the end of inflation when the spacetime can be considered 
a vacuum (which has an upper bound of total number of e-foldings for inflation). This 
power spectrum can be seen as an approximation to the curvature perturbation power 
spectrum, and the vacuum here can be seen as the vacuum with respect to curvature 
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perturbations. Length scales of interest for the the CMB anisotropics give N k of order 
of 60, and therefore one expects the theoretical ambiguity on the power spectrum to be 
sizeable if the total duration of the de Sitter stage corresponds to a number of e-foldings 
not far from 60. Of course, in this case it would be difficult to assume that the spacetime 
is in a vacuum state (see for example [19, 20, 21]). 

3.2 Adiabatic vacuum 

In situations in which the vacuum is defined through the notion of particles, the natural 
fundamental operator is the number operator which counts the number of particles with 
momentum k. Parker in Ref. [11] postulated certain reasonable conditions that the 
number operator for a scalar field must satisfy in a FRW universe. The conditions were 
as follows: 

1. N k be Hermitian due to the counting interpretation. 

2. When the expansion is stopped at any time (i.e. a/a — 0), the operator becomes 
the usual Minkowski number operator. 

3. The vacuum expectation value of the number operator varies slowly as possible 
with time as the expansion rate a/ a becomes arbitrarily slow. 

The first two conditions are obviously reasonable, and give rise to the definition of the 
number operator as 

N k {n) ee a™ a™ (20) 

where the superscripted r x refers to the time at which the boundary condition for vacuum 
is set (or equivalently, the boundary conditions for hk{r)). 

The third condition is the statement that the vacuum should be defined as to keep 
the number of particles as unchanging as possible. This third condition is what Parker 
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has called the minimization postulate [11] and was later developed further by Parker and 
Fulling [12]. The formalism developed based on satisfying these conditions is called the 
adiabatic vacuum formalism [12, 8]. 

It is important to note that to define the number operator N k , one must define what 
one means by a particle at any given time, T\. On the other hand, to have defined a 
vacuum at an earlier time t < T\ (recall that Heisenberg representation states such 
as the vacuum are time independent), one needed a definition of a particle at time r . 
Hence, the present formulation requires that one define particles at two different times. 
This should be contrasted with the formalism of computing the vacuum expectation 
value of the stress tensor, where the definition of particles (or equivalently the vacuum) 
needs to be defined only once, since the stress energy tensor is a local quantity which 
does not rely on the basis of Fourier expansion. However, since we will be taking the 
same prescription for particles at time To and time T\ in such a way that the total number 
of conditions that has to be specified is the same as in the case of computing the vacuum 
expectation value of the stress tensor, this inherently nonlocal definition of vacuum will 
probably also minimize the growth of the stress tensor to a large extent. Nonetheless, it 
is not obvious how the conclusions of our analyses would differ if stress tensor vacuum 
expectation values are used instead of N k in carrying out the minimization postulate. 
We will defer this question to a future work and focus on the particle based formalisms 
in this paper. 

The adiabatic formalism specifies the value and the first time derivative of the mode 
function tit at a fixed time for any fixed momentum. The boundary condition data is 
specified by matching to an asymptotic expansion (adiabatic expansion) of the mode 
equation to all orders in the asymptotic parameter (even though the expansion does not 
converge in general, by the very nature of an asymptotic expansion). In de Sitter space, 
for example, the infinite adiabatic order vacuum is identical (up to the uncertainties 
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inherent to the adiabatic vacuum) to the Bunch-Davies vacuum. In inflationary space- 
times in which the scale factor cannot be made arbitrarily small due to the finite duration 
of the inflationary phase or in which the momentum cannot be made arbitrarily large 
due to a cutoff, the adiabatic formalism naively has a chance to still specify a unique 
vacuum, unlike the Bunch Davies prescription. However, we will find that the adiabatic 
prescription suffers from an ambiguity problem as well. 

The construction of the adiabatic vacuum is as follows. First, define the concept 
of an adiabatic order as the power of 1/T that results for any term in an 1/T — > 
asymptotic expansion after one makes the transformation r — > r and d/dr — > T~ l d/dr 
in the differential equation for the modes h k . 7 ■ Then, any prototypical mode equation of 
the form 

Mr) + ^(r)Mr)=0 (21) 

with w\ = k 2 + m 2 a 2 {r) + (6£ — l)a/a (where £ is a constant and m is the mass of (f>) 
turns into 

l^t k (r)+w 2 h k (r) = (22) 

where wl = k 2 + m 2 a 2 (r) + (6£ — l)d/(T 2 a) and a tilde has been added to h k to be a 
reminder that this function carries a fictitious parameter T (later this will be set to 1 at 
which time the function will be denoted as hk)- Now, make a change in variables from 
hk to Wk by writing 

(where the T in the exponent should be noted) and from Eq. (22), we obtain a new 
differential equation 

W, 3 /li/,A 2 

(24) 



W 2 k =w\- 



2T 2 



fwA 2 

W k 2\W k 



7 Note that the point about only the derivative being transformed has been missed by a recent paper 
[23]. 
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We can then define a map 



A 



Wl 



\ 



Wk 2T 2 



wf 



3 / WP 



2 \ 



(25) 



which is a map that raises the adiabatic order by two and also define 



n+2] 



= A 



wl 



(26) 



where the superscript denotes the adiabatic order and = a//c 2 + m 2 a 2 . All of this 
construction lead to an approximate mode equation solution good to Ath adiabatic order 
in asymptotic expansion in 1/T — > of 



{^fK A \r')dr'T) 



2W, 



[A] 



(27) 



This asymptotic expansion solution will serve as a template, just to set the boundary 
conditions for the mode function hk which can in general be written as 



h k (r) = A k f k (r) + B k f* k (r) 



(28) 



where f k are exact basis functions satisfying Eq. (22). Specifically we define the Ath 
adiabatic (order) vacuum at time r by using the boundary condition 

Kir,) = h [ k A] (r ) + O(l/T^), 



(29) 



where the left-hand side is the exact mode solution to the prototypical differential 
equation Eq. (22) and the boundary conditions are enforced to only order \/T^ A+ ^ 
as T — > oo. 8 After the boundary conditions are set, one sets the vacuum mode of Eq. (2) 

as 



K{t) = K{t)\ t=1 



(30) 



8 The correction is order l/r( j4 + 1 ) and not of order 1/T( A + 2 ) because there is an integration in the 
exponent. 
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which removes the fictitious parameter T. 

To see that this construction satisfies the minimization postulate if we take the vac- 
uum adiabatic order A to oo, let us compute the vacuum expectation value of the num- 
ber operator corresponding to the number density per mode. This requires a Bogoliubov 
transformation from the vacuum mode solution with the boundary condition at r = tq 
into the one with the boundary condition at r = T\ (any later time at which the particles 
are no longer being created). Defining the Bogoliubov transformation as 

h^(r) = a k hl (T)+(3 k h^(r) (31) 

the number density of particles per momentum k is given by n k = \Pkf- If the vacuum 
in the past is defined at r = t with infinite adiabatic order boundary condition and the 
vacuum today is defined at r = T\ with infinite adiabatic order boundary condition, car- 
rying out the Bogoliubov transformation with the solution written in the form Eq. (23), 
one finds 

1 f 1 / W T ° W T1 \ 2 1 
ro (0|^( Tl )|0) To oc l/^ro)! 2 = |- - + (W? - W-) 2 j(32) 

where the right hand side (composed of exact solutions to Eq. (24)) can be evaluated 
at any r and the superscripts indicate the time at which the boundary conditions were 
placed. 

Now, the minimization postulate is satisfied if J^r(r o (0|A' r fc(ri)|0) ro ) is minimized for 
all non-negative integers n if the expansion rate can be turned off arbitrarily slowly. For 
any a(r), we can affect the slowly turning off of the expansion rate by reintroducing the 
adiabatic order parameter l/T — > as before. Since by construction, the WJ) functions 
in Eq. (32) when expanded in l/T match the asymptotic expansions in l/T of 
of Eq. (26), as long as the asymptotic expansion is uniform in r between r and T\ (no 
singularities in the asymptotic expansion occur between To and ti), the right hand side 
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of Eq. (32) vanishes identically when expanded in 1/T. Hence ^r(r O (0|-/V fc (ri)|0) r() ) will 
fall off faster than any finite power of 1/T as T — > oo, thus satisfying the minimization 
condition. 

In general, one must remember that does not exist because the construc- 

tion procedure generates an asymptotic expansion about a nonanalytic point rather than 
a convergent series. In other words, for a fixed order A (and with T = 1), there exists 
only an A-dependent region in time for which approximates well the exact solution 
to Eq. (5), with the leading error on the approximation growing with A for any fixed 
time r. Hence, when the limit A — > oo is taken first, the time region in which the 
approximation is valid can shrink to for a fixed T — 1. 

As an example of an adiabatic vacuum, let us apply this formalism to the case of 
massless scalar field in dS space (considered in Eq. (4)). The iteration map Eq. (25) 
produces up to sixth order 



W, 



[0] 2 



= k 2 



Wl 2]2 = k 2 



[4] 2 



k 2 



Wf 2 = k 2 



1 - 
1 - 
1 



(Tkrf 
2 



+ 



{Tkr) 2 (Tkr) 4 
2 3 



+ o 



(TkrY 



O 



(33) 



(Tkr) 2 (Tkr) 4 (T/cr) 6 \(Tkr) 

The template asymptotic expansion h\f^ obtained from approximates the exact 

solution with an error of order 1/T A+1 . For example, one can write down explicitly 
using wjf^ (dropping the higher order corrections to it) 

-ikrT r A 2l 



h [4] 



1 



krT + 5£; 5 r 5 T 5 



+ 0(1/T b 



(34) 
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where the 1/T 5 term is the leading uncertain term (displayed just for clarity), and up to 
1/T 4 , the template function matches the first term of the exact solution 

-ikrT / ■ \ ikrT / • \ 

h ^ A ^{ l -WT)) +B >wA 1+ ^)- (35) 

Carrying out the matching procedure of Eq. (29) except up to 5th adiabatic order instead 
of the just the required 4th, one finds 

* - 1+ 2^ + - (36) 
B k = + C(l/T 6 ) (37) 

where we have displayed one of the uncertain terms explicitly for clarity. Note that the 
coefficients in general depend on r (when the boundary condition was placed), but the 
sensitivity to it is to higher adiabatic order. In general, the coefficients A k and B k will 
be of the form 

A k = 1 + (9(1/T A+1 ), 

B k = + O{l/T A+1 ) (38) 

for an Ath order adiabatic vacuum. In the limit that A — > oo, the uncertainty drops off 
faster than any finite power of 1/T. 

Recall that in the Bunch-Davies prescription, one could not remove the uncertainty 
in the choice of {A k , B k } if we did not have information to let a(r) — > 0, say because of 
the finite period of inflation. Since in the adiabatic prescription, the boundary conditions 
may be set to infinite adiabatic order at a finite time, let us see whether at finite initial 
time T it is possible to define an infinite adiabatic vacuum without any ambiguity. By 
inspection, we can cast 2 in the form 



A/2 / 1 \" / 1 

^]=» = ^2(-l)»(„ + l)^_j + o< 



(39) 
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This means that in the limit A — > oo, 2 is expressed as 



W W = k 2 ^:(-ir(n + l)(-±-) (40) 



n=0 



KkrT 



2n 



which converges to a simple function 9 

W [oo] 2 = p {krTf 

(l + (fcrT) 2 ) 2 1 J 

Miraculously, the series converges in this simple situation! 10 Indeed, inserting into 
(27), we obtain at any time r 

= 9r ( x " ) ' (42) 

which in view of Eq. (29) will set A k = 1 and = in Eq. (35) (note that we have 
normalized the function properly with the arbitrary integration constant freedom). 

This seems to indicate that the vacuum has been uniquely specified. Unfortunately, 
this is not true since if an infinite adiabatic order vacuum with T — > oo is chosen to be 
hk(r), we can always choose another vacuum 



h k = yfl + \B k \*h k {r) + B k h* k (r) (43) 

if B k h* k (T) falls off faster than any finite power of 1/T. The infinite order adiabatic 
vacuum boundary conditions do not distinguish h k and h k . Indeed, this can be seen 
directly in Eq. (29) because the "equation" is asymptotic, up to terms that vanish faster 
than Ath power of 1/T, where in the infinite adiabatic order case A = oo: i.e. we should 
have written instead of Eq. (42), the equation 

—ikrT / ' \ 

ht ] = b^[ 1 -k?T) + ° {1/TOO) (44) 



9 Notice that in the limit r = t — > — oo one recovers W^°°J 2 = k 2 , i.e. the infinite adiabatic order 
vacuum reduces to the BD vacuum. 

10 We will later give an example where there is no convergence. 
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where 0(1/T°°) indicates the nonperturbative term possibly dropped in the matching 
to the template function. Since T follows every factor of k, we can guess the uncertainty 
by assuming an exponential function for B k as 

B k » exp(-fcT/(a(ro)^(r ))) (45) 

which should be compared to Eq. (12). (Note that this would also serve as a good 
estimate of uncertainties even outside of inflationary phase.) This leads to 



~ C(exp(-e JV °- 7Vfe )), (46) 



V^k) 

where we have used the same notation as Eq. (19). Because of the double exponential, 
the uncertainty is very quickly negligible as N becomes larger than N k . 

In the case in which the expansion in 1/T stops converging, as in the quasi-dS case, 
we may be able to estimate the uncertainty in the adiabatic vacuum formalism a little 
less arbitrarily (compared to Eq. (45)) as follows. Indeed, if at some adiabatic order 
n* the expansion in 1/T stops converging (with T = 1, the higher order terms give a 
larger correction to the mode solution than the n*th order term), this signifies that the 
nonadiabaticity in the system enters at the n*th derivative. Therefore, as far as assigning 
an uncertainty to the adiabatic vacuum is concerned, instead of using Eq. (45), we may 
estimate the uncertainty in fixing the vacuum to simply be the term in the asymptotic 
expansion. Although the uncertainty obtained this way is not "nonperturbative" , it is not 
clear what the physical advantage is in artificially "turning off" the expansion (using 1/T) 
beyond the n* derivative which is in some sense the limit of adiabaticity characteristic 
of the physical system. We will later use this approach to compute the uncertainty for 
the slow roll inflationary case. 

There is an exception to the existence of an uncertainty in the adiabatic vacuum 
formalism, which is already suggested by the estimate in Eq. (45) in the case that a(r ) = 
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O. 11 Namely, if the boundary conditions are set at a time in which the nonadiabaticity 
is identically 0, then there is no expansion in 1/T, and hence there is no ambiguity 
in the adiabatic vacuum. Indeed, the loss of time translational invariance which is 
operationally at the heart of the uncertainty of the vacuum disappears at the time when 
the nonadiabaticity is identically 0. The test of the disappearance of the nonadiabaticity 
is that for any A becomes identical, or more explicitly 



for all A where To is the time at which the vacuum is defined. This time To typically is ±00, 
which means that an asymptotic expansion in r must be taken to match the boundary 
conditions. Hence, the reason why the ambiguity disappears in this exceptional case can 
also be seen as due to the asymptotic expansion being in time r instead of 1/T, since 
the coefficients in Eq. (28) are time independent although they can be be T dependent. 

This exceptional case can be viewed as a meeting point of the adiabatic vacuum 
formalism and the Bunch-Davies formalism. Its spirit is similar to the Bunch-Davies 
formalism in that the asymptotic expansion can be in time r instead of a fictitious 
adiabatic parameter T. However instead of matching on to a Minkowski prescription, it 
matches on to a zeroth adiabatic order WKB prescription as in the adiabatic vacuum 
formalism. Essentially, the adiabatic formalism states that although not in general, 
sometimes a unique vacuum can be defined in regions of spacetime where the frequency 
is adiabatically a constant to an arbitrarily good degree just as in the Bunch-Davies 
prescription, except with the frequency given by the leading order WKB ansatz instead 
of a Minkowski prescription (Minkowski prescription means that with the scale factor 

frozen at a fixed value). 

11 Here, we are not claiming a (to) = gcnerically defines an adiabatic point, but merely that Eq. (45) 
suggests there are special points in spacetime where the uncertainty vanishes. 




(47) 
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The reader should be aware that the adiabatic formalism is compatible with the 
Bunch-Davies formalism. For example, this formalism always gives the same vacuum as 
the Bunch-Davies formalism whenever the Bunch-Davies formalism can define a unique 
vacuum. Also, we would like to comment that although we would like to extend the 
"uniqueness" to the case in which WKB iteration Eq. (25) converges for n = oo (i.e. 
W^°°\r) is well defined with T fixed at a finite value), it is not obvious to us whether there 
is any reason to consider such situations in any special manner since the convergence of 
WKB iteration is not identical to the absence of nonadiabaticity. It merely means that 
the asymptotic expansion in 1/T is being taken about an analytic point in the function 
generating the asymptotic expansion. 

4 An Example of Calculable Nonperturbative Par- 
ticle Production 

In this section, we will discuss a model from Ref. [24] which illustrates the adiabatic 
formalism in situations where nonperturbative particle production can be unambiguously 
calculated. Owing to the supposed background metric, the vacuum construction will be 
unique, and hence exponentially suppressed particle production will be calculable. 

Before moving on to the discussion, we would like to alert the reader here that this 
model is also reviewed on pg. 70 of Ref. [8]. However, their usage of adiabatic prescription 
can be misleading (it was to us, at least). As a consequence of their loose usage, they 
give the wrong impression that they can calculate reliably an exponentially suppressed 
particle production even without taking to ±oo the time at which the vacuum is specified. 
For example, one sees that the matching in their equation (3.117) involves taking A — > oo, 
which by itself results in a possibly exponentially suppressed ambiguity in the vacuum 
that we discussed in Eq. (45). However, as long as the vacuum is specified at r = — oo, 
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we agree with their results modulo their inconsequential phase error in their equation 
(3.124). We present our reanalysis here such that the reader would not be confused 
about the discrepancy of how Ref. [8] seems to claim that nonperturbative calculations 
can be unambiguously done in the context of the adiabatic vacuum formalism even 
without taking tq — > — oo while we claim that such calculations cannot be done, strictly 
speaking, without being able to take such a limit for the time at which the vacuum is 
specified. 

To see how the nonperturbative Bogoliubov coefficient arises in an exactly solvable 
model of Ref. [24], consider the conformally coupled mode equation (i.e. Eq. (21) with 
£ = 1/6) where the scale factor is given by 

a 2 (r)=a 2 + (3 2 T 2 (48) 

leading to the effective frequency 

w\ = ml3\ + m 2 l3 2 r 2 (49) 

where 

ma 2 k 2 

x = — + ^p- (50) 

Since we would like to carry out the adiabatic procedure, we scale the derivatives w.r.t. 
r by 1/T: 

(^+^)^ = 0. (51) 

For differential equations of the form 

^ + (x 2 + A) M = 0, (52) 
dx z 

the general solution is the parabolic cylinder function 

= D - { i«wW + i>] (53) 
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u 



and its complex conjugate where for our case, one should identify 

x = ^m(3Tr , (54) 

A = TA , (55) 

where one notes that as T — > oo, the solution is effectively doubly scaled as \/Tt and 
TA. Note also that the parabolic cylinder function D p (z) is defined by its boundary 
condition 

D p (z) ~ z p e~ z2/4 (56) 

as z — > oo. However, this asymptotic expansion is not the same expansion as the T — > oo 
expansion since there is an effective double scaling of r and A. 

Let us first see what the adiabatic vacuum procedure gives for the particle production. 
We begin by writing explicitly the exact mode function (solution to Eq. (51)) as 

h k = A k f + B k f* (57) 

where 

7.1/4 

m = j^^ e ' nXT/8D -(^)/^ - (58) 

which as one sees has the curious double scaling of r and A mentioned before. We 
need to obtain an asymptotic expansion of the parabolic cylinder function D p (z) as 
T — > 00, which corresponds to an asymptotic expansion as the complex order p goes 
to ioo. Unfortunately, we do not know of any systematic way of doing this because 
asymptotic expansion for D p known in the literature corresponds to either the limit 

\jml3T\r\ > AT/2 (59) 



or 



AT/2 > m(3TT 2 , (60) 
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both of which are not satisfactory for generic r. However, for r — ► ±00 we can satisfy 
Eq. (59), and for r^Owe can satisfy Eq. (60). In reproducing the result of Audretsch 
and Schaeffer, we will focus on the limit r — > ±00 satisfying Eq. (59). 

In the limit of r — > —00, the parabolic cylinder function has the asymptotic formula 



D p (z) ~ z p e~ z2/A 



where 



P 



x r 



1 - ]T 0(p 2n /z 2n ) + 

n=l 

2 ~l/4+i\T/4 e n\T/8+in/8 ( m ^) ^ (1-iAT) 

' i-(i-iAT)/2 ™£Zr 2 



"e 2 



(AT) 



2n 



71=1 



{m(3T) n T 



n~-2n 



z = (i — 1)J mf3Tr , 



;i-iAT)/2, 



(61) 



(62) 



(63) 
(64) 



giving 



/ 



e 2 



=2 aT/ V 7r/8 (m/3T) iAT/4 1 r \ iXT/2 



where one has to be careful to write 



(65) 



(i - l)y/m/3TT = -{i - l)^Jml3T\r\ (66) 
for r < 0. 12 

Now, we construct the template function. Let us start with the zeroth order expansion 
Wf\r) = \/ m p\ + m 2 f5 2 r 2 



which gives 



h 



[0] 



[2(v^r + VA + /?mr a )]- iAT / 2 



exp 



-iT 



/3m (A + f3mr 2 )T 



(67) 



(68) 



12 We see that in Eq. (62), we cannot take T — > 00 after taking r — > —00 as the T expansion is out of 
control in that limit (i.e. T limit and r limit do not commute unless the T expansion can be summed 
up). 
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In the limit r — > — oo, we therefore have 



2m/5|r| 



exp 



.m/3r|r|T 



—2- 



(^3|r|)^/ 2 (A)- 



■iXT/2 



(69) 



We can consider higher orders: 



Wl 2] = J(3m{\ + /3mr 2 ) + 



(3/3mr 2 - 2A) v //?m(A + /?mr 2 ) 



8T 2 (A + /3mr 2 ) s 



[i] 



/9m(A + /9mr 2 ) + 



(3/?mr 2 - 2A) v //3m(A + (3mr 2 ) 
8T 2 (A + /5mr 2 ) 3 



(76A 2 - 732/3 Amr 2 + 297/5 2 m 2 r 4 ) J/5m(A + /5mr 2 ) 



128T 4 (A + /3mr 2 ) 6 
In the limit that r — > — oo, these have the expansions 



(70) 

(71) 
(72) 



/3to|t| + 



AM 



A 2 |r| 



2r 2 8/3mr z 



+ ... +7^ 



3|t| 



+ 



ir , , 1 ( ~297|r| 



+ 



19A|r| 

T 2 V8/9mr 4 ' 16/? 2 m 2 r 6 
4731AM \ 



+ 



(73) 
(74) 



T 4 \ v 128/5 3 m 3 r 8 256 / 3 4 m 4 r 10 
where we note that each power of 1/T n has an expansion in while the leading 

correction coming from 1/T n terms is \r\/ \\f3m] a ~ l T 2n ). This means that even for a 
fixed finite T, we can write 



lim W ia] {t )/W [0] {t ) 

Tq— > — OO 



(75) 



and hence we can conclude that this corresponds to an exceptional situation in which 
infinite adiabaticity occurs at the time t = — oo when the vacuum is set. At any finite 
time To, we would have to use the matching in expansion of 1/T, which can lead to 
nonperturbative ambiguity, but here, all the terms in 1/T are effectively identically 0. 
For completeness, we integrate to display the template function exponent. We find 



T J" dxWl 4] {x) = 



1 „ m ATlrl , 
-(3Tm\r\r + ' ' In \r\ + ... 

2 2r 



1 

+ T 



t \Qf3mT 2 



+ ... 



+ o (t- 3 ) 



(76) 
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which means that when exponentiated, both the 1/T and 1/r expansions are both well 
in control. 

Because of Eq. (75), not only do Eqs. (69) and (65) match up to r independent phases 
(which can always be chosen appropriately in the indefinite integral Eq. (27)), but the 
specification is unique without any nonperturbative ambiguity. Furthermore, we have 
from Eq. (75) that h k °\r ~ — oo) = h k °°\r ~ — oo), implying that the matched vacuum 
h k is an infinite adiabatic order vacuum. We shall denote this vacuum as 



(77) 



where we have denoted the placement of the boundary condition at r = — oo. 

Since we can read off the template function at r = oo from Eqs. (68) and (69), we 
can write 



r) = hr-°°(-ry 



(78) 



specifying the vacuum at a future time. By examining the asymptotic expansions 9.246 
of [25], we can then write 



h?=-°°i 



giving 13 a particle production coefficient of 



r > 0) 



h T k 1= °°i 



r) + 



-re 



-ttA/2 



h T k 1= °°* 



(r), 



(79) 



fa 



(80) 



• -^TA/2 



More explicitly, if we had not set T = 1 through Eq. (30) we could write f3 k = —ie 

which shows explicitly that this coefficient is to all orders in 1/T. Therefore, this 

nonperturbatively suppressed particle production could not have been reliably calculated 

if we had to make an expansion in 1/T at any point in the calculation (we have made 

an expansion in 1/T in this section just to show that it was not necessary as can be seen 
13 There is an error in the book by Birrcll and Davies [8] for the first coefficient. 
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for example in Eq. (75)). Such an expansion in 1/T would have been necessary if the 
time r at which vacuum is specified were finite. 



5 The adiabatic vacuum in the slow roll spacetime 

For any realistic scenario of inflation, the background spacetime is quasi-dS, unlike the 
exact dS that we treated explicitly before. By quasi-dS, we mean the Hubble rate 
(H = d(r)/a 2 (r) ) is not exactly constant, but changes with time as H/a = —eH 2 , where 
e is one of the so-called slow-roll parameters. 14 The other useful slow-roll parameter is 
V = (V" /3H 2 ) where V(4>) is the inflaton potential and primes denote derivatives with 
respect to the inflaton field 0. We wish to estimate the uncertainty in the vacuum in 
such a spacetime, and thereby estimate the uncertainty in the inflationary prediction of 
density perturbations. 

In the gauge invariant treatment of cosmological perturbations generated during in- 
flation, it is useful to define the gauge invariant Mukhanov variable u = aScj)^ + 
[26] where z = <j)/H and 5(f)(°^ \1/ are the gauge invariant inflaton fluctuation and the 
Bardeen gravitational potential respectively [27]. Performing for the variable u a field 
expansion similar to Eq. (1), one finds the mode equation (after inserting the adiabatic 
parameter 1/T) 



h k (r) = 0, (81) 



T 2 r 2 

where p = 3(3e — i]). The exact mode equation solution can be written in terms of Hankel 
functions as 

h k = ^V^T [A k H^\-kTr)e^ + ^ + B k H^\-kTr)e-^ + ^] (82) 

where one must be careful to note that — r > (we have also taken out Tkr independent 
14 Recall we are using conformal time coordinates. 
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phases from Ak and for convenience and normalized conveniently). Here, the order 
is defined as v = y / 9~+4p/2. 15 

Let us now construct the adiabatic vacuum following the recipe of subsection 3.2. 
After rescaling d/dr — > {l/T)d/dr in the mode equation, the template frequencies can 
be constructed as 



Wl 



Wl 2] 2 (r) 



w, 



[4] 2, 



k 



r) = k 2 



1 - 
1 - 



(2+p) 
(krT) 2 

(2+p) 



(fcrT) 



+ c 



2 + 1 3+ 2^(^tF + ° 

1 - (2+p) (^F + ( 3 + ^)(^ 

' 1 5 " 



(krTY 



(83) 
(84) 

(85) 



+ -4 



TP 



(A;tT) c 



(ArrT)* 



(86) 



However, unlike in the dS case, the series is generally non-convergent for any finite r = t . 
As we discussed before, this type of nonconvergence is more generic than the convergent 
case. 

To fourth order in adiabatic expansion, we find the following WKB template function 

+ p/2) p ip 5p \ 



-ikrT 



1 + 



'2k \ kTr Ak 2 T 2 r 2 6k 3 T 3 r 3 24A; 4 T 4 r 4 / 

For the asymptotic expansion of the Hankel functions H^' 2 \ we have from Ref. [25] 



(87) 





^ (-l)* r(i/ + fc + |) (-i)»r(z/ + n + | 

^ (2^) fe fc!r(i/ — A; + |) 1 (2i*) B r (i/ - n + |) 



1 I> + A; + ±) „ 1 T(v + n + \ 



,(88) 



E 

k=0 



(2iz) k k\T{v - k 



(2w) B r (i/ - n + i) 



,(89) 



15 The Hankel's functions are defined as = J v + iY v and = *, where J v and Y„ are, 
respectively, the Bessel functions of the first and the second kind. 
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where Re(V) > —1/2, |arg(z)| < n and 9 li2 are coefficients smaller than unity in front of 
the remainders. A useful formula is 

I> + A; + ±) _ (4z/ 2 - l 2 )(4z/ 2 - 3 2 )...(4z/ 2 - (2k - l) 2 ) 



Y(v-k + \) 2 2fe 
Expanding to first order in p and fourth order in T, one can then compute 



(90) 



h k (r) « 



(a*/ + B k n + ^{-^ (i + f) / + 3 fe (i + f ) ./" ; 



A(krT) 2 " hJ kJ J QikrTf 

(91) 



+ 24^Tf {Akf + Bkf * } 



where / = e~ tkrT . Matching ^ and h$ just by inspection (comparing Eqs. (87) and 
(91) ), one sees that A k — 1 and 5 fc = as expected. 

Just to check further, suppose we took the 4th order template h$ and expanded to 
6th order in T" 1 and matched it to h k to 6th order (instead of just to 4th order). To 
accomplish this, note that since only needs to be expanded to T~ 6 accuracy even 
though it is being multiplied by T since the next order term in expansion is T~ 8 , 
resulting in a correction of T~ 7 which we are ignoring anyway. We find 

-ikrT r _ 3i 9 + 5p 



^(6th) -h k + ^j- 
This results in 



8A; 5 r 5 T 5( 4 + 3 ^ + 4A;6 r 6 T 6 



(92) 



Ak = 1 + %f^ + ^ 6 ^ (93) 



15(2+p) 2 . fcTro 



where r is the time at which the boundary conditions are placed. We see explicitly 
the order of the residual corrections to the 4th order asymptotic expansion, and it is as 
expected. One can easily check similarly that the 6th adiabatic order WKB template 
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function also results in A k = 1 and B k = to 6th order in 1/T. More explicitly, we have 
the template function 

a? 1 = 4 41 + 



(95) 



r Sip 7p 



Wt 5 T 5 8A; 6 t 6 T 6 

that fixes the vacuum. Indeed, the fact A k — 1 and T4 = is no surprise since the for- 
mulae Eqs. (88) and (89) and the template functions are the same asymptotic expansions 
to any given order. 

The matching at time r and at any given adiabatic order A of the template function 
h k A ^ with the mode function (82) suffers of an ambiguity introduced by the nonvanishing 
remainders (proportional to 9i) displayed in expansions (88) and (89). As discussed ear- 
lier, the order at which the remainder stops converging (with T = 1) can be attributed 
to be the extent to which the vacuum is nonadiabatic. Hence, instead of a guessed non- 
perturbative uncertainty of Eq. (45), we can set T = 1 and compute the uncertainty that 
is due to the intrinsic nonadiabaticity of the background spacetime. (This is something 
we discussed at the end of Sec. 3.2.) 

More explicitly, we can look for the adiabatic order n* at which the (absolute values 
of the) remainders are minimized, and then say that we cannot define the vacuum more 
certainly than this order due to the nonadiabaticity. Expanding the remainders to first- 
order in the slow-roll parameters, we can express the absolute value of the remainder 
as 



R{n) = 



-l) n r(z/ + n+ 1/2) 



(2ik\r\) n r(u -n + 1/2) 



pT(2 + n)T(n- 1) 



+ 0(p 2 ) . (96) 



3 • 2»(fc|7b|)" 

We see clearly that the remainder vanishes in the dS limit of p — 0, which corresponds 
to e = rj — 0. This is another way of understanding why in the pure de Sitter case we 
have been able to construct a convergent W^ 00 '. At a formal level, it merely means that 
the asymptotic expansion point 1/T — > is an analytic point, allowing a convergent 
asymptotic expansion. 
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We can find the minimum of R at the location by taking a discrete derivative i?(n* + 
1) — R(n*) = 0. Although will not lie at an integer value in general, the actual solution 
will be at the nearest integer. 16 We thus find 



integer (-^ + 1^1 + 8(1 + fc|r |)) « integer ^/2fc|r |) (97) 



to be the order at which the uncertainty is minimized where the integer function finds 
the nearest integer value of the argument. Note that the value is independent of p to 
leading order because in the limit that p vanishes, the corrections identically vanish, 
giving an "arbitrary" n*. In terms of e-folds, we can write 



y/2exp((N -N k )/2) (98) 



where Nk denotes the number of e-foldings from the time when a given wavelength 
A = a/k leaves the horizon (A = 1/H) till the end of inflation and N denotes the 
number of e-foldings from the time the vacuum was set to the end of inflation. As we 
already mentioned, length scales of interest for the the CMB anisotropies give N k of 
order of 60. 

The corresponding ambiguity in the value of the parameter can be calculated from 
R(n) as follows. We can consider the template function to be 



h [ r 1] = \ m.HW(-kT)e^i+ rn ., (99) 



where r n is related to the remainder as 



r n = hfR^n), (100) 

ikr 

"° = 73T (101) 
= vs ( ir r( :+;i + i)e-^> 

V ; {2i\r\k) n r(u-n + \) V ; 



16 When solving, however, one will find it convenient to assume self-consistently that n* will in the 
end be at an integer value. 
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where |-Ri(n)| = \/2R(n). Now, solving for the coefficients A k and B k using the system 

4"*~Vo) = ^(to)|t=i, (103) 
h [ r 1] '(ro) = Kfo)\T=i (104) 

gives 

B k « -ftH 2 ^ . (105) 
ar 

The ambiguity in the power spectrum Vn of the comoving curvature perturbation 1Z = 
u/z is given by 



2 \ReB k \ p exp(-2 v / 2e (iVo - JVfc)/2 - (iV - N k )/2) , (106) 



where we are using the notation and approximation explained below Eq. (19). Note that 
this is generically slightly larger than Eq. (46) and that it artificially vanishes in the limit 
that p — > 0. In the case that p — » 0, the best estimate that we can give for the uncertainty 
is Eq. (46), which is a merely a guessed function that drops off nonperturbatively fast. 
The theoretical ambiguity on the power spectrum of the comoving curvature perturbation 
is sizeable if the total duration of the de Sitter stage corresponds to a number of e-foldings 
not far from 60. Of course, in this regime, it may not be a good approximation to treat 
the spacetime to be void of fluctuations. 

6 Conclusions 

Within the context of Bunch-Davies vacuum formalism and the adiabatic vacuum for- 
malism we have answered the following question: "What is the minimal theoretical 
uncertainty in the inflationary curvature perturbation calculation if we assume that the 
curvature perturbation state at sometime r near the beginning of inflation was a slow 
roll vacuum?" Even without any trans-Planckian effects, effective field theory cutoff 
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related effects, or nongravitational field interaction effects, there is a minimal uncer- 
tainty in the curvature perturbation predictions of inflation coming from the inability 
to uniquely specify a vacuum (due to gravitational interactions). Within the adiabatic 
vacuum formalism, the power spectrum uncertainty is given by Eq. (106), and in more 
general situations (applicable to outside of the inflationary regime) by Eq. (45). The 
Bunch-Davies formalism gives a larger uncertainty of Eq. (19). 

The minimum uncertainty presented here applies to all of the efforts [13, 14, 21, 22] 
to obtain measurable small effects on the CMB. In practice of comparing with data, 
there will certainly be other theoretical uncertainties, not only from other interaction 
effects of the inflaton, but reheating historical uncertainties [28] as well as astrophysical 
uncertainties which will most likely overwhelm the minimal uncertainty presented here, 
unless the number of e-foldings is very close to the minimum required for inflation (e.g. 
around 60). Even in that case, however, most likely, it will be difficult to assume that the 
curvature perturbation quantum state is that of a vacuum due to other energy density 
fluctuations present which depend on the history leading to the initiation of inflation. 
Nonetheless, it is important and reassuring to know that the inflationary vacuum in the 
conservative sense has very little ambiguity contrary to the impressions given particularly 
by Ref. [14]. Indeed, the main qualitative conclusion one can draw from our work is 
that even with a finite period of inflation, we can in most cases neglect the uncertainty 
associated with the vacuum when defined according to the most reasonable particle based 
definitions. 

Finally, we would like to comment that the general idea that the vacuum is uncertain 
with a finite period of inflation is not entirely new (see for example [7, 29]). Here, we 
merely tried to quantify this in the context of slow roll inflationary models and show 
explicitly that the best estimates for the uncertainty lead us not to worry about this 
effect unless the total duration of inflation is close to 60 e-folds or so. 
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